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Abstract 

Many neutrino mixing scenarios that have \x — r symmetry with #13 = are 
in disagreement with recent experimental results that indicate a nonzero value for 
#13. We investigate the effect of small perturbations on Majorana mass matrices 
with n — t symmetry and derive analytic formulae for the corrections to the mix- 
ing angles. We find that since m\ and m2 are nearly degenerate, fj, — r symmetry 
mixing scenarios are able to explain the experimental data with about the same 
size perturbation for most values of #12. This suggests that the underlying unper- 
turbed mixing need not have #12 close to the experimentally preferred value. One 
consequence of this is that a new class of models with fj, — r symmetry is possible, 
with unperturbed #12 equal to zero or 90° for arbitrary unperturbed $13. 
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Of the numerous neutrino mixing scenarios discussed in the literature [T] , several have 
\l — t symmetry, such as tri-bimaximal mixing (TBM) [2], bimaximal mixing (BM) [3J, 
hexagonal mixing (HM) [4] and scenarios of A§ mixing [5]. In these scenarios, #23 = 45°, 
#13 = 0, and only #12 depends on the particular model. Tri-bimaximal mixing is most 
popular because the value of 6>i 2 predicted by TBM is close to that preferred by the 
current experimental data. However, the latest results from the T2K [6], MINOS [7], 
and Double Chooz [8J experiments suggest a nonzero value of #13, and the recent Daya 
Bay [9] and RENO [10J experiments find 9 13 ^ at the 5.2<r and 4.9<r level, respectively. 
Various corrections may reconcile such models with nonzero #13 pQ. In this Letter we 
consider small perturbations acting on Majorana mass matrices with /i — r symmetry 
and estimate the size of perturbations required to explain the experimental data. 

We find that for /i — r symmetries with almost any initial value of # 12 (i.e., before 
the perturbation), the minimal size of the perturbations needed to bring the model in 
agreement with experimental data varies by only about 20%. The reason is that the #12 
correction depends only on the ratio of perturbation terms and not on their absolute 
size, and the overall size of the perturbation is determined by the corrections to #13 and 
#23, which are relatively small. We also show that a new category of models with fi — r 
symmetry, #23 = 45°, 612 = or 90°, and arbitrary #13, can also fit the data with small 
perturbations. 

We start with the mass matrix for Majorana neutrinos 

M = U*M di£Lg rf , (1) 

where M diag = diag(mi, m2, m 3 ), U is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
mixing matrix [TT] (without the multiplicative diagonal matrix of Majorana phases), and 
we work in the basis in which the charged lepton mass matrix is diagonal. The masses 
m2 and 777.3 are complex and mi can be taken to be real and non-negative. 

The general condition describing [l — r symmetry (also sometimes called fi — r uni- 
versality) is [12j 

\U^\ = \U Ti \, for 1 = 1,2,3. (2) 
From the standard form of the mixing matrix these conditions are equivalent to 

9 23 = 45° , Re(cos 6 12 sin 12 sin 6i 3 e iS ) = . (3) 

Hence, there are four classes of // — r symmetry: (a) #23 = 45°, #13 = 0; (b) #23 = 
45°, #12 = 0; (c) 9 23 = 45°, 9 12 = 90°; (d) 9 23 = 45°, 5 = ±90°. Class (a) contains models 
with tri-bimaximal, bimaximal, hexagonal, and symmetries, while class (d) includes 
tetramaximal symmetry [13] . Classes (b) and (c) have not been studied before because 
the unperturbed #12 angle is far from the experimentally preferred value, but, as we show 
below, small perturbations can have a large effect on #12, and therefore these models 
should not be ignored. 
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Class (a): 0° 3 = 45° , 6° n = 



We first examine the effect of small perturbations on models in class (a). The initial 
(unperturbed) mixing matrix can be written as 



/cosfl?. 



Un = 



12 



sin ^ 2 



J 12 



V2 
sinfl° 2 



V2 

cos 6»J 2 



1 

^2 



(4) 



\/2 V2/ 



and the initial mass matrix is 
M = U*Mf^Ul = 



m\c\ 2 + m°s? 2 
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^2 
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where MQ iag = diag(m ( j ) , m 2 , m"), and Cjfc, Sjfe denotes cos0j k and sin 0° fe respectively. 
Under a small perturbation the final (resultant) mass matrix can be written as 



M = U*M* iag U j +E, 
where the perturbation matrix E has the general form 



E = M-M, 



o 



e ll e 12 ^13 
^12 ^22 ^23 
v Cl3 ^23 ^33, 



(6) 



(7) 



Treating the three masses as eigenvalues of the mass matrix with each column of 
the mixing matrix as the corresponding eigenvector, we can use traditional perturbation 
methods to find the corrections to the three angles and three masses. From experiment 
we know that mi and rri2 are nearly degenerate, so that degenerate perturbation theory 
with | (Sm^ | ^ l^ m 3il an d < \ m< l\ (where Sm^ = mj — m", and the index k denotes 
the heaviest eigenstate), can be used. For simplicity, we assume Mq and E are real and 
employ the following notation: 



e l — e ll j ^2 — €i2 + Ci3 , €3 — €12 — Ci3 , 64 — €22 + ^33 + 2f23 j 
^5 — e 22 — ^33 , ^6 = £22 + £33 — 2f23 — 2en . 

We find the first order corrections to the three masses to be 



(8) 



4ei+e 6 ± (^25m° 21 -^8e 2 3 + e| + A(5m° 21 ) 2 + 45m§ 1 (2v / 2e 3 sm29° 12 + e 6 cos2fl° 2 ! 



where the plus sign is for i = 1 and the minus sign is for i = 2, and 
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(i) 



1 



-e 4 . 



(9) 
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Table 1: Best-fit values and 2a ranges of the oscillation parameters [H] used to find the 



€ij, with 6m 2 = 



m 2 \ 



— ml and Am 2 = 



m 3 



| 2 - (ml + 



m 2 \ 2 )/2. 



Hierarchy 


M°) 


M°) 


M°) 


r5m 2 (10- 5 eV^) 


|Am 2 |(l(T 3 eV 2 ) 


Normal 


33.6+20 


8.91JS 


38.4+^ 


7 

' - d ^-0.39 


2 43+ u - i2 


Inverted 


33.6+££ 




38.81^ © 47.5 - 53.2 


' - J ^-0.39 


2 42 +aii 



The first order corrections to the mixing angles are 



60S 



69% 



69 
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and the second order correction to 9i 2 is 



69 



(2) 
12 



V2e 2 e 5 cos 2{6° 12 + 69%) + (e 2 - e 2 /2) sin 2 
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Imposing l^ml^l <C l^m^-J, the expressions for an d ^li/ simplify to 



3(1) 



66$ ~ 



65 



2<5mjj 1 ' 



V2e2_ 
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(15) 



We note that while 69% and (Wjg' 1 are suppressed by a factor of order tj/6m® l , to leading 
order 5# 12 depends only on ratios of linear combinations of e%, e§ and 6m 2l (which is 
0{eij)). Therefore large corrections to 9\ 2 are possible even for small corrections to #23 
and 6> 13 . 

A recent global three-neutrino fit [14] yields the parameter values in Table [TJ We 
have done a numerical search to find perturbed mass matrices that give the oscillation 
parameters and which have small perturbations. In our search, we first fix #23 = 45° and 
9® 3 = 0, consistent with /i — r symmetry, and choose a particular value for 9® 2 and the 
magnitude of m\ for the normal hierarchy (or m 3 for the inverted hierarchy). The global 
fit in Table [TJ then defines the magnitudes of the other two final masses and the three 
final mixing angles (since 9\ 3 = 0, the initial Dirac phase does not matter). 

We characterize the size of the perturbation as the root-mean-square (RMS) value of 
the perturbations, i.e., 



3(1) 
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Table 2: Top half: values of the perturbations (in 10~ 3 eV) that give the best-fit param- 
eters in Tabled] and have the minimum €rms f° r the given 9® 2 , for the normal hierarchy 
and mi = 0. Bottom half: representative values that fit the experimental data within 
2a and for which all have a similar magnitude (with mj = 0, m 2 = 0.0054 eV, 
ml = 0.0595 eV, m 1 = 0.0072 eV, 5 = 180° and all other phases equal to 0). 





en 


ei2 


da 


^22 


e 2 3 


^33 


£rms 


60 


-3.05 


-3.50 


-5.99 


-2.72 


-1.52 


5.77 


4.10 


45 (BM) 


-1.32 


-4.74 


-4.74 


-3.58 


-0.66 


4.90 


3.79 


35.3 (TBM) 


0.32 


-4.66 


-4.82 


-4.40 


0.16 


4.08 


3.74 


30 (HM) 


1.07 


-4.31 


-5.18 


-4.78 


0.54 


3.71 


3.79 





0.00 


-1.38 


-8.11 


-4.24 


0.00 


4.24 


4.36 


60 


5.41 


-4.17 


-4.52 


-5.00 


-9.94 


3.36 


6.14 


45 (BM) 


6.76 


-4.43 


-4.26 


-5.67 


-9.27 


2.69 


6.08 


35.3 (TBM) 


7.66 


-4.32 


-4.37 


-6.12 


-8.82 


2.24 


6.08 


30 (HM) 


8.11 


-4.17 


-4.52 


-6.35 


-8.59 


2.01 


6.09 





9.46 


-2.52 


-6.17 


-7.02 


-7.92 


1.34 


6.28 



where i and j sum over neutrino flavors. Hence, srms is determined by the following 
quantities: three initial masses, two initial Majorana phases, two final Majorana phases 
and one final Dirac phase. We scan over these quantities with all phases taken to be 
either or 180° to find the minimum value of erms for a given 9® 2 . We follow the same 
procedure for classes (b) and (c) below. For class (d), all values of the phases are allowed. 

We show the perturbations that give the smallest €rms fo r the normal hierarchy, 
mi = and several values of m Table El It is clear that the sizes of brms a re 
approximately the same regardless of the value of 9® 2 ; we find that the smallest €rms 
for each 9® 2 varies by at most 17% for the examples shown. This can be explained 
by the perturbation results derived above as follows. From Eq. flH]) we have £rms = 

\J e\ + e\ + e| + |e§ + \e\ + \{2e\ + e§) 2 /3; since m 3 ^> m 1 ,m 2 for the normal hierarchy 
with mi = eV and the first order perturbations of the three masses are much smaller 
than m.3, we can assume Sm^ ~ m\ ~ m% m V Am 2 = 0.0493 eV. Then from Eq. ( I15jl 
we know that in order to get the correction <5#23 = —6.6° and <5#i3 = 8.9° for any value of 
0? 2 , we need e 5 = -0.0114 eV and e 2 = 0.0108 eV, so that y/ef+ 4/2/3 = 0.00449 eV, 
which is already close to the e^MS values found in Table [2j The small discrepancy can 
be explained by the perturbation of the three masses and other e's. Hence, we can say 
that the size of the perturbation mainly comes from the corrections to #23 and #13. From 
Eq. ( fill we know that the correction to 9 u is determined by the relative ratio of 63 to t§ 
and the actual size of the perturbation does not matter. This means that we can have 
large corrections for 9\ 2 with a (relatively) small perturbation. 

We note that initial values of 9 12 on the "dark side" {9\ 2 > 45° and m\ < m 2 ) can 
also fit the data with perturbations that are similar in magnitude to those needed for 
tri-bimaximal mixing (see the entry for 9® 2 — 60° in Table EJ). 
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Table 3: Top half: same as Table |2j except for the inverted hierarchy and m 3 = 0. 
Bottom half: same as Table [21 except for the inverted hierarchy and m\ = 0.05 eV, 
m" = 0.052 eV, m° = 0, m 3 = 0.002 eV, and all phases equal to 0. 





en 


ei2 


ei3 


^22 


£23 


e 33 


CRMS 


60 


-0.86 


-4.94 


-5.64 


5.57 


-0.43 


-4.72 


4.31 


45 (BM) 


-0.47 


-5.29 


-5.29 


5.38 


-0.23 


-4.91 


4.29 


35.3 (TBM) 


-0.05 


-5.30 


-5.28 


5.17 


0.03 


-5.12 


4.28 


30 (HM) 


0.16 


-5.23 


-5.36 


5.07 


0.08 


-5.22 


4.28 





0.00 


-4.47 


-6.12 


5.15 


0.00 


-5.15 


4.32 


60 


-3.56 


-4.89 


-5.27 


5.95 


2.67 


-3.92 


4.49 


45 (BM) 


-3.06 


-4.98 


-5.18 


5.70 


2.92 


-4.17 


4.47 


35.3 (TBM) 


-2.73 


-4.94 


-5.22 


5.54 


3.08 


-4.34 


4.46 


30 (HM) 


-2.56 


-4.89 


-5.27 


5.45 


3.17 


-4.42 


4.46 





-2.06 


-4.28 


-5.88 


5.20 


3.42 


-4.67 


4.50 



In the top half of Table [21 en and €23 are much smaller than the other for some 
values of 9\ 2 . We have checked that if these values are set to zero, the experimental 
constraints can still be satisfied at the 2a level without a large change in the nonzero 
parameters. Therefore if some perturbations are exactly zero due to symmetries, the 
resulting mass matrix can still fit the experimental data with small perturbations. 

For the inverted hierarchy, some representative sets of that give the minimum crms 
are shown in Table [3] for 777.3 — 0. The minimum e^MS as a function of 9\ 2 varies only by 
about 1% in this case, i.e., the minimum €rms varies with e° 12 even less for the inverted 
hierarchy than for the normal hierarchy. 

Clearly, if perturbations are large enough that tri-bimaximal mixing can explain the 
experimental data, then other // — r mixing scenarios, such as bimaximal, hexagonal 
mixing and A 5 mixing, can also explain the experimental data with about the same size 
perturbation. Hence, tri-bimaximal mixing has no special position among the fi — r 
symmetry mixing scenarios when a perturbation is required to fit the experimental data. 
Also, it is possible for all the perturbations to have a similar magnitude and still give the 
oscillation parameters within their 2a ranges; see the bottom half of Tables [2] and [3j 

We also varied the size of the final masses by changing the value of mi in the normal 
hierarchy and m 3 in the inverted hierarchy. We find that the minimum e^MS decreases 
as the size of the final masses increases for both the normal and inverted hierarchies. 
For the quasi-degenerate hierarchy (in which the magnitude of the absolute masses is 
larger than V Am 2 ) the size of the perturbation can be very small. This can be explained 
by the perturbation equations: since <5m 31 ~ m 3 — mi Am 2 /(m 3 + mi) for small 
perturbations, and Am 2 is fixed by experimental data, then 5m 31 will decrease if the 
masses increase, and similarly for 5m 32 . Then Eqs. ( IT21) and ( IT31) show that in order to 
get the same corrections for and 8 23 , the size of the perturbation should also decrease. 
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Classes (b) and (c): 8 2: 



23 



45°, 6° 12 = or 90° 



For class (b) (#23 = 45°, 9\ 2 = 0); since the Dirac phase is irrelevant, the initial mixing 
matrix and mass matrix can be written as 



/cos ^ 3 



U 



V2 
_ sin e° 3 

\ 7T 
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V2 
1 
V2 



sm 1 



13 



cos 0j 3 
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cos 6>?„ 
V2 I 



(17) 



and the initial mass matrix is 



M 



/ 777,9 



0„2 



m l c 13 + m s s 



S^IS 



V 



(mjj-?ng)si3Ci3 
72 



rm 



+ m °3 C 13 



+ m l S l3J 



\{-ml 



m° 3 c 2 13 + m?s? 3 ) 



(mg-mg)si 3 ci 3 
V2 



\ 



(18) 



If we redefine the phase of the wavefunction 7^3 to — ^3, or change the initial angle #23 
from 45° to 135° and switch the indices 2 and 3, then the mass matrix in Eq. (fl8|) is 
exactly the same as that in Eq. 

For the above initial mass matrix, corrections must shift 812 from to 33.6°, and # 13 
from the initial arbitrary angle to 9.0°. We used the same scan procedure as before and 
searched for the minimum €rms f° r various values of 8® 3 (see Table H]). We find that 
for 6® 3 < 20°, the data can be explained with about the same size perturbation as was 
found for class (a). For example, when 6® 3 = for class (b), the initial mass matrix is 
the same as 9i 2 = for class (a), and therefore the minimum srms is also the same. 
In particular, when 8® 3 is close to 9.0° in class (b), the minimum €rms is even smaller 
than the minimum value for class (a) because the correction to #13 is smaller in this case. 
Although the correction to 6*12 is large, it does not affect the size of the perturbation too 
much because its size is mainly determined by the corrections to #13 and #23, as noted 
before. However, for 59i 3 greater than about 20°, the size of the perturbation required 
to fit the data becomes larger since #13 must change by more than 10°. 

For class (c) (#23 = 45°, 9\ 2 — 90°), we find that switching with m 2 makes the 
initial mass matrix the same as the initial mass matrix of class (b). Since we scan all 
possible values of m? and m 2 , the minimum €rms for a given 6® 3 for class (c) is the same 
as for class (b). 



Class (d): 6% 3 = 45°, 5° = ±90° 



If we fix e° 23 



45°, 5° = ±90° and vary both 9® 2 and ^^3, this category includes mixing 
scenarios such as the tetramaximal mixing pattern (T 4 M) [T3] , and the correlative mixing 
pattern with 5 = ±90° [15]. For 6»° 3 < 20° and 6° 12 < 45°, the smallest e RMS for the normal 
hierarchy (with mi = 0) varies from 2.29 x 1(T 3 eV to 5.26 x 10" 3 eV, where the minimum 



value occurs at 9l 3 



9° and & 



12 



32°, and the maximum value occurs at 



'13 



20° and 
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Table 4: Top half: same as Tabled except for class (b) (9\ 2 — 0). Bottom half: same as 
Table [21 except for class (b). 





en 


ei2 


ei3 


e22 


e23 


e33 


£rms 





0.00 


-1.38 


-8.11 


-4.24 


0.00 


4.24 


4.36 


5 


0.48 


1.44 


-5.28 


-4.48 


-0.24 


4.00 


3.27 


10 


-0.44 


4.21 


-2.52 


-4.02 


0.22 


4.46 


3.06 


15 


-2.64 


6.59 


-0.14 


-2.92 


1.32 


5.56 


3.90 


20 


-5.85 


8.30 


1.57 


-1.32 


2.93 


7.17 


5.24 





9.46 


-2.52 


-6.17 


-7.02 


-7.92 


1.34 


6.28 


5 


9.01 


1.13 


-2.52 


-6.80 


-7.69 


1.56 


5.41 


10 


7.66 


4.67 


1.02 


-6.12 


-7.02 


2.24 


5.22 


15 


5.47 


8.00 


4.35 


-5.03 


-5.93 


3.33 


5.80 


20 


2.50 


11.00 


7.35 


-3.54 


-4.44 


4.82 


6.92 



6® 2 = 0. Therefore small perturbations can fit the experimental data for a wide range of 
9i 2 an d 6*i3 for class (d). 

In summary, we studied small perturbations to Majorana mass matrices with fi — r 
symmetry that yield experimentally preferred oscillation parameters. We find that the 
size of the perturbations (which decreases as the neutrino mass scale is increased), is 
mainly determined by the corrections to #23 and #13, and that small perturbations can 
give a very large correction to #12 because to first order, the #12 correction depends only 
on the ratio of perturbation terms and not on their absolute size. Hence, most mixing 
scenarios with fi — r symmetry can explain the experimental data with perturbations of 
similar magnitude, and tri-bimaximal mixing has no special place among scenarios with 
fi — t symmetry. We also find that slightly perturbed \i — r symmetric models with 
12 = or 90° are viable for 13 < 20°. 
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